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ABSTRACT 

We report on the analysis of the Kepler short-cadence (SC) hght curve of V344 Lyr obtained during 
2009 June 20 through 2010 Mar 19 (Q2-Q4). The system is an SU UMa star showing dwarf nova 
outbursts and superout bursts, and promises to be a touchstone for CV studies for the foreseeable 
future. The system displays both positive and negative superhumps with periods of 2.20 and 2.06-hr, 
respectively, and we identify an orbital period of 2.11-hr. The positive superhumps have a maximum 
amplitude of ^0.25-mag, the negative superhumps a maximum amplitude of '^0.8 mag, and the orbital 
period at quiescence has an amplitude of ^0.025 mag. The quality of the Kepler data is such that 
we can test vigorously the models for accretion disk dynamics that have been emerging in the past 
several years. The SC data for V344 Lyr are consistent with the model that two physical sources 
yield positive superhumps: early in the superoutburst, the superhump signal is generated by viscous 
dissipation within the periodically flexing disk, but late in the superoutburst, the signal is generated 
as the accretion stream bright spot sweeps around the rim of the non-axisymmetric disk. The disk 
superhumps are roughly anti-phased with the stream/late superhumps. The V344 Lyr data also reveal 
negative superhumps arising from accretion onto a tilted disk precessing in the retrograde direction, 
and suggest that negative superhumps may appear during the decline of DN outbursts. The period 
of negative superhumps has a positive P in between outbursts. 

Subject headings: novae, cataclysmic variables; dwarf novae; stars: individual (V344 Lyr); white 
dwarfs; hydrodynamics 



1. INTRODUCTION 

Cataclysmic variable (CV) binary systems typically 
consist of low-mass main sequence stars that transfer 
mass though the LI inner Lagrange point and onto a 
white dwarf primary via an accretion disk. Within the 
disk, viscosity acts to transport angular momentum out- 
ward in radius, allowing mass to move inward and ac- 
crete onto the primary white d warf (e.g. iWarnedFiggSal : 
iFrank et al.ll2002t lHelliedl200l"l ). In the case of steady- 
state accretion the disk is the brightest component of 
the system, with a disk luminosity Ldisk ~ GMiMi/ Ri, 
where Mi is the mass accretion rate onto a white dwarf 
of mass Ml and radius Ri . 

While members of the novalike (NL) CV subclass dis- 
play a nearly constant mean system luminosity, members 
of the dwarf nova (DN) subclass display quasi-periodic 
outbursts of a few magnitudes thought to arise from a 
thermal instability in the disk. Specifically, models sug- 
gest a heating wave rapidly transitions the disk to a hot, 
high- viscosity state which significantly enhances Mi for 



|wood@fit.edu' 

^ Department of Physics and Space Sciences, Florida Institute 
of Technology, Melbourne, FL 32901, USA 

2 NASA Ames Research Center, Moffett Field, CA 94095, 
USA 

^ Bay Area Environmental Research Institute, Inc., 560 Third 
St. West, Sonoma, CA 95476, USA 

* National Optical Astronomy Observatory, Tucson, AZ 
85719, USA 

^ CRESST and Astroparticle Physics Laboratory 
NASA/GSFC, Greenbelt, MD 20771, USA 

® Department of Physics, University of Maryland, Baltimore 
County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 

^ NASA/Goddard Space FUght Center, Greenbelt, MD 20771, 
USA 



a few days. Furthermore, within the DN subclass there 
are the SU UMa systems that in addition to normal DN 
outbursts display superoutbursts which are up to a mag- 
nitude brighter and last a few times longer than the DN 
outbursts. The SU UMa stars are characterized by the 
appearance at superoutburst of periodic large- amplitude 
photometric signals (termed positive superhumps) with 
periods a few percent longer than the system orbital pe- 
riods. So-called negative superhumps (with periods a few 
percent shorter than Porb) are also observed in some SU 
UMa systems. 

The oscillation modes (i.e., eigenfrequencies) of any 
physical object are a direct function of the structure of 
that object, and thus an intensive study of SU UMa su- 
perhumps that can make use of both a nearly-ideal time- 
series data set as well as detailed three-dimensional high- 
resolution numerical models has the potential to eventu- 
ally unlock many of the long-standing puzzles in accre- 
tion disk physics. For example, a fundamental question 
in astrophysical hydrodynamics is the nature of viscos- 
ity in differentially rotating plasma disks. It is typically 
thought to result from the magnctorotational instabil- 
ity (M RI) proposed bv ,Balbus fc H awlcv (1998 . ): BalbJ^ 
(|2003f l. but the observations to-date have been insuffi- 
cient to test the model. 

1.1. V344 Lyrae 

The Kepler field of view includes 12 CVs in the Ke- 
pler Input Catalog (KIC) that have published results 
at the time of this writing. Ten (10) o f thes e sys- 
tems are listed in Table 1 of iStill et al.l ()2010l . here- 
after Paper I). Two additional systems have been an- 
nounced since that publication, the dwarf nova sys- 
tem BOKS-45906 (KIC 9778689) (,Feldmeier et al...20111 
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and the A M CVn star SPSS J1 90817.07+394036.4 (KIC 
4547333) (iFontaine et alJ[20TTI) . 

The star V344 Lyr (KIC 7 65957 ] is a SU UMa star 
that Hes in the Kepler field. iKatol ()1993l ) observed the 
star during a superoutburst {V ~ 14), and reported the 
detection of superhumps w ith a period P = 2.1948 ± 
0.0005 hr. In a later study iKato et all H20Q2D reported 
that the DN outbursts have a recurrence timescale of 
16 ± 3 d, and that the superout b ursts have a recurrence 
timescale of '^110 d. lAk et all (|2008f ) estimated a dis- 
tance of 619 pc for the star using a period-luminosity 
relationship. 

In Paper I we reported preliminary findings for V344 
Lyr based on the second-quarter (Q2) Kepler observa- 
tions, during which Kepler observed the star with a '^1- 
min cadence, obtaining over 123,000 photometric mea- 
surements. In that paper we reported on a periodic sig- 
nal at quiescence that was either the orbital or negative 
superhump period, and the fact that the positive super- 
hump signal persisted into quiescence and through the 
follo wing dwarf nova outburst. 

In iCannizzo et al.l (|2010L hereafter Paper II) we pre- 
sented time-dependent modeling based on the accretion 
disk limit cycle model for the 270 d (Q2-Q4) light curve 
of V344 Lyr. We reported that the main decay of the 
superoutbursts is nearly perfectly exponential, decaying 
at a rate of '^12 d mag~^, and that the normal outbursts 
display a decay rate that is faster-than-exponential. In 
addition, we noted that the two superoutbursts are ini- 
tiated by a normal outburst. Using the standard accre- 
tion disk limit cycle model, we were able to reproduce 
the main features of the outburst light curve of V344 
Lyr. W e significantly expand on this in ICannizzo et al.l 
(j2011| ) where we present the 1-year outburst properties 
of both V344 Lyr and V1504 Cyg (Cannizzo et al. 2011). 

In this work, we report in detail on the results ob- 
tained by studying the Kepler Q2-Q4 data, which com- 
prise without question the single-best data set obtained 
to-date from a cataclysmic variable star. The data set 
reveals signals from the orbital period as well as from 
positive and negative superhumps. 

2. REVIEW OF SUPERHUMPS AND EXAMPLES 

Before digging into the data, we briefly review the 
physical processes that lead to the photometric modu- 
lations termed superhumps. 

2.1. Positive superhumps and the two-source model 

The accretion disk of a typical dwarf nova CV that is 
in quiescence has a low disk viscosity and so inefficient 
exchange of angular momentum. As a result, the mass 
transfer rate Mli through the inner Lagrange point LI is 
higher than the mass transfer rate Mi onto the primary. 
Thus, mass accumulates in the disk until a critical surface 
density is reached at some annulus, and the fl uid in that 
annu l us transitions to a hi gh-viscosity state (|Cannizzol 
Il998t ICannizzo et al.ll2010D . This high- viscosity state 
propagates inward and/or outward in radius until the en- 
tire disk is in a high- viscosity state characterized by very 
efficient angular momentum an d mass transport - the 
standard DN outburst (see, e.g.. lCannizzolll993atlLasotal 
120011 for reviews). In this state. Mi > Mli and the disk 
drains mass onto the primary white dwarf. 



During each DN outburst, however, the angular mo- 
mentum transport acts to expand the outer disk radius 
slightly, and after a few to several of these, an otherwise 
normal DN outburst can expand the outer radius of the 
disk to the inner Lindblad resonance (near the 3:1 coro- 
tation resonance). This can onl y occur for systeni s with 
mass ratios q = M2/M1 < 0.35 (jWood et al.il2009[) . 

Once sufficient mass is present at the resonance radius, 
the common superhump oscillation mode can be driven 
to amplitudes that yield photometric oscillations. The 
superhump oscillation has a period P+ which is a few 
percent longer than the orbital period, where the frac- 
tional period excess e+ is defined as 
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These are the are the so-called common or positive su- 
perhumps, where the latter term reflects the sign of 
the period excess e+. In addition to the SU UMa 
stars, positive superhumps have also been observed in 
novalike CVs ^ Patterson et al."1993'; Retter et al."1997t 
i Skillman et al.l J997; Patterson ct al. 2005; Kim et al] 

l2009f ). the interacting bin ary wh ite dwa r f AM 

CVn sta r s (iPatterson et all |l993t l\yarneH Il995bl: 
iNelemansI 120051 : iRoelofs et al.l 120071: IFontaine et al ' 



I2OIII), and in low-mass X-ra y binaries fCharles et al 
1991; M ineshige ct al 119921: FO 'Donoghuc & Charles 
ll996tlRetter et al.ll200a iHvnes et al.il20060 . 

Figure [T] shows snapshots from one full orbit of 
a smoothed particle hydrodynamics (SPH) simulation 
{q = 0.25, 100,000 particles) as well as the associ- 
ated simulation light curve ('seelSimpson fc WoodI [19981: 
IWood fc Bur"yi2007l : IWood et al.ll2009t) . The disk parti- 
cles are color-coded by the change in internal energy over 
the previous timestep, and the Roche lobes and positions 
of Ml are also shown. Panels 1 and 6 of Figure [T] shows 
the geometry of the disk at superhump maximum. Note 
that here the superhump light source is viscous dissipa- 
tion resulting from the compression of the disk opposite 
the secondary star. The local density and shear in this 
region are both high, leading to enhanced viscous dissi- 
pation in the strongly convergent flows. The orbit sam- 
pled in the Figure is characteristic of early superhumps 
where the disk oscillation mode is saturated, and the re- 
sulting amplitude significantly higher (^^0.15 mag) than 
the models produce when dynamical equilibrium (~0.03 
mag) due to the lower mean energy production in the 
models at superhump onset. 

As a further detail, we note that whereas the 2 spiral 
dissipation waves are stationary in the co-rotating frame 
before the onset of the superhump oscillation, once the 
oscillation begins, the spiral arms advance in the pro- 
grade direction by ~180° in the co-rotating frame dur- 
ing each superhump cycle. This prograde advancement 
can be seen by careful inspection of the panels in Figure 
[TJ Indeed, this motion of the spiral dissipation waves 
is central to the superhump oscillation - a spiral arm 
is "cast" outward as it rotates through the tidal field 
of the secondary, and then brightens shortly afterward 
as it compresses back into the disk i n a converging fiow 
(|Smith et al.ll2007l: IWood et al.ll2009l) . 

While viscous dissipation within the periodically- 
flexing disk provides the dominant source of the super- 




Fig. 1. — Snapshots from one orbit of a g = 0.25 SPH accretion disk simulation shortly after superhump onset. The 100,000 simulation 
particles are color-coded by their internal energy change ("luminosity") over the previous simulation time step. The simulation light curve 
shown is calculated as the sum of the luminosities of all the particles. Note that superhump maximum corresponds to frames 1 and 6, 
where the disk opposite the secondary is radially compressed and hence a region of strongly convergent flows. Note also that the accretion 
stream impact region is deeper in the potential well in frame 3 than in frames 1 and 6 - this yields the "late" superhump signal. 



hump modulation, the accretion stream bright spot also 
provides a periodic photometric signal whe n sweeping 
around the rim of a non -axisy mmetric disk ()Vogtl[l982t 
[Osaki 1985t iWhitebl^ 1988; Kunzell2004[ ). The bright 
spot will be most luminous when it impacts most deeply 
in the potential well of the primary (e.g., panel 3 of Fig- 
ure [U and fainter when it impacts the rim further from 
the white dwarf primary (Panels 1 and 6). This signal 
is swamped by the superhumps generated by the flexing 
disk early in the superoutburst, but dominates once the 
disk is significantly drained of matter and returns to low 
state. The disk will continue to oscillate although the 
driving is much diminished, and thus the stream mecha- 
nism will continue to yield a periodic photometric signal 
of decreasing amplitude until the oscillations cease com- 
pletely. 

This photometric signal is what is termed late su- 



perhumps in th e liter a ture (e . ; 



Patterson et"all [2000l l2002al: J l 

Sterken et al. f l2007tlKato eral1[2009. ' .201Ql 



iHessman et al.l 119921; 

Tcmplcton et al. 200a 

.. iRolfe et ahl 

()2001|) presented a detailed study of the deeply eclipsing 
dwarf nova lY UMa observed during the late superhump 
phase where they fo und exactly this beh avior. They used 
the shadow method iWood et al.l (|1986( ) to determine the 
radial location of the bright spot (disk edge) in 22 eclipses 
observed using time-series photometry. They found that 
the disk was elliptical and precessing slowly at the beat 
frequency of the orbital and superhump frequencies, and 
that the brightness of the stream-disk impact region var- 
ied as the square of the rel ative velocity of the stream 
and disk material (see also iSmakI l201Cl() . Put another 
way, the bright spot was brighter when it was located on 
the periastron quadrant of the elliptical disk, and fainter 
on the apastron quadrant. 
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Thus, two distinct physical mechanisms give rise to 
positive superhumps: viscous dissipation in the flexing 
disk, driven by the resonance with the tidal field of the 
secondary, and the time-variable viscous dissipation of 
the bright spot as it sweeps around the rim of a non- 
axisymmetric disl^. For the remainder of this paper we 
refer to this as the two-source m odel of positive super- 
humps (see also lKunzil2002l |2004[) . These two signals are 
approximately antiphased, and in systems where both 
operate at roughly equal amplitude, the Fourier trans- 
form of the light curve can show a larger amplitude for 
the second harmonic (first overtone) than for the funda- 
mental (first harmonic). 

As an example of this double humped light curve, in 
Figure[2]we show 20 orbits of the q = 0.25 simulation dis- 
cussed above (Figure [T]) starting at orbit 400, by which 
time the system had settled into a state of dynamical 
equilibrium. The inset in this Figure shows the the av- 
erage superhump pulse shape obtained from orbits 400- 
500 of the simulation, where we have set phase zero to 
primary minimum. Note that here the average pulse 
shape is complex but approximately double-peaked. The 
Fourier transform displays maximum power at twice the 
fundamental frequency. When we examine the disk pro- 
files, we find that the dominant peak arises from the 
disk superhump described above, but the secondary peak 
roughly half a cycle later results from the impact of the 
bright spot deeper in the potential well of the primary 
(see panel 4 of Figure 1). The substructure of this sec- 
ondary maximum results from the interaction of the ac- 
cretion stream with the spiral arm structures that ad- 
vance progradely in the co-rotating frame. Panel 3 of 
Figure [T] is representative of the disk structure at the 
time of the the small dip in brightness observed at su- 
perhump phase 0.55. The dip is explained by the fact 
that the accretion stream bright spot at this phase is lo- 
cated in the low-density inter-arm region, and therefore 
that the accretion stream can dissipate its energy over a 
longer distance. In addition the oscillating disk geometry 
results in this region having a larg er radius, and lower 
velocity contrast near this phase. iHowell et all ()1996f ) 
discuss the observation and phase evolution of the two 
secondary humps in the SU UMa system TV Corvi. 

The 3 AM CVn (helium CV) s ystems that are in per - 
manent high state - AM CVn (jSkillman et al.l 119990 . 
HP Lib (jPatterson et a l.' '2002a') and the system SDSS 
J190817.0 7-h394036.4 (KIC 004 547333) announced re- 
cently by iFontaine et al.l ()201ll ) - all display average 
pulse shapes that are strongly double humped. AM 
CVn itself is frequently observed to show no power in 
the Fourier transform at the fundamental superhump 
oscillation frequency (|Smak 1967; Faulkner et al. 1972; 
[Patterson et al.l 119921: ISkillman et al.l 119991) . AM CVn 
systems are known to be helium mass transfer systems 
with orbital pe riods ranging between 5 min and ^1 hr 
(see reviews bv IWarned[l995bl: ISolheiml[2010l) . 

In contrast, the hydrogen-rich old-novae and novalike 

* For completeness, we note that recently Smak (2009, 2011) 
has proposed that the standard model, described above, does not 
explain the physical source of observed superhump oscillations. In- 
stead, he suggests that irradiation on the face of the secondary 
is modulated, which yields a modulated mass transfer rate A/li, 
which in turn results in modulated dissipation of the kinetic energy 
of the stream. 
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Fig. 2. — Twenty orbits of the simulation light curve once the 
dynamical equilibrium state has been reached in the g = 0.25 sim- 
ulation discussed above. The light curve is double humped. The 
inset shows the average light curve calculated from orbits 400 to 
500, with the phase set to zero for primary minimum. Primary 
maximum again corresponds to modulation in the disk as discussed. 
Secondary maximum is caused by the accretion stream impact spot 
hitting minimum radius of the non-axisymmetric disk. The sub- 
structure of the secondary maximum is discussed in the text. 

CVs that show permanent superhumps display mean 
pulse shapes that are nearly always similar to the satu- 
ration phase light curves as shown in Figure [TJ and there 
is no example we know of where a permanent superhump 
system shows a strong double-humped light curve. The 
reason for this is clear upon reflection: the AM CVn 
disks are physically much smaller than the disks in sys- 
tems with hydrogen-rich secondary stars, resulting in a 
much higher specific kinetic energy to be dissipated at 
the bright spot since the disk rim is much deeper in the 
potential well of the primary. The smaller disk may also 
yield a smaller amplitude for the disk oscillation signal. 
In the hydrogen-rich systems in permanent outburst, the 
disks are large, the mass transfer rates are high, and the 
disk signal dominates, with a relatively minor contribu- 
tion from the stream source. 

We tested the viability of the two-source model 
through three additional numerical experiments. First, 
we again restarted the above simulation at orbit 400, 
but now with the accretion flow through LI shut off 
completely. In this run, there is no accretion stream 
and hence no bright spot contribution. We show the 
first 20 orbits of the simulation light curve in Figure [3l 
With the stream present, the light curve has the double- 
humped shape of Figure[2]above, but without the stream 
the light curve is sharply peaked with no hint of a dou- 
ble hump. Note that because there is no low-specific- 
angular-momentum material accreting at the edge of the 
disk, the disk can expand further into the driving zone. 
This expansion results in the pulse shape growing in am- 
plitude as the mean disk luminosity drops. The pulse 
shape averaged over orbits 410-440 is shown as an inset 
in the Figure, and clearly shows that the oscillating disk 
is the only source of modulation in the light curve - max- 
imum brightness corresponds to a disk geometry like that 
from panel 1 of Figure [T] above. The mean brightness is 
roughly constant for orbits 410-440, and at orbit 440 the 
mean brightness and pulse amplitude begin to decline as 
some 50% of the initially-present SPH disk particles are 
accreted by orbit 450. 

Our second test was to restart the simulation a third 
time at orbit 400, but this time to enhance the injec- 
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Fig. 3. — The first 20 orbits after SPH simulation restart with no 
mass flow through LI. The pulse shape increases in amplitude as 
the mean luminosity increases, resulting from the expansion of the 
disk. Note that the pulse shape is strongly - and singly - peaked. 
The mean pulse shape calculated over orbits 410-440 is shown in 
the inset. 
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Fig. 4. — The flrst 20 orbits after SPH simulation restart with 
enhanced mass flow through LI as described in the text. The 
pulse shape is saw-toothed, and maximum light corresponds to 
the phase where the bright spot is deepest in the potential well of 
the non-axisymmetric flexing disk. The average pulse shape over 
these 20 orbits is shown as an inset. The mean system luminosity 
increased with the simulation particle count. The mean pulse shape 
calculated over orbits 410-440 is shown in the inset. 

tion rate of SPH particles (mass flow) at LI by roughly 
a factor of 2 over that required to keep the disk par- 
ticle count constant (Figure S]). This enhanced mass 
flux again dramatically changes the character of the light 
curve. Here the mean pulse shape as shown in the inset 
is saw-toothed, but with the substructure near the peak 
from the interaction of the stream with the periodic mo- 
tion of the spiral features in the disk as viewed in the 
co-rotating frame. Careful comparison of the times of 
maximum in these two runs (Figures [3] and |4]) reveals 
that they are antiphased with each other. For example, 
the simulation light curve in Figure [3] shows maxima at 
times of 403.0 and 404.0 orbits, whereas the simulation 
light curve in Figure |4] shows minima at these same times. 

Our third experiment was more crude, but still ef- 
fective. We began with a disk from a g = 0.2 low- 
viscosity SPH simulation run that was in a stable, non- 
oscillating state. We offset all of the the SPH par- 
ticles an amount 0.03a along the line of centers [i.e., 
{x,y,z) {x + O.G3a, y, z)], scaled the SPH particle 
speeds (but not directions) using the vis viva equation 



(2) 



Fig. 5. — A snapshot at orbit 401.35 from a. q = 0.40 SPH sim- 
ulation showing negative super humps. The green line running di- 
agonally though the primary indicates the location of the line of 
nodes; the disk midplane includes this line, but is below the orbital 
plane to the right of the line, and above the orbital plane to the 
left of the line. The particles are color coded by luminosity, and 
the brightest are rendered larger. It can be seen that the accretion 
stream impact region at this phase is deep in the potential well of 
the primary, roughly at the line of nodes. The deeper the accretion 
stream impact region, the brighter the bright spot. The simulation 
light curve is derived from the "surface" particles calculated using 
a simply ray-trace technique. 

and restarted the simulation. This technique gives us 
disk which is non-axisymmetric but not undergoing the 
superhump oscillation. The results were as expected: we 
find maxima in the simulation light curves at the phases 
where the accretion stream impacts the disk edge deepest 
in the potential well of the primary. 

In summary, numerical simulations reproduce the two- 
source model for positive superhumps. 

2.2. Negative Superhumps 

Photometric signals with periods a few percent shorter 
than Porb have also been observed in several DN, nova- 
likes, and AM CVn systems - in some cases simulta- 
neously with posi tive superhumps (see, e.g.. Table 2 of 
iWood et al.|[200l and Woudt et al. 2009). These os- 
cillations have been termed negative superhumps owing 
to the sign of the period "excess" obtained using Equa- 
tion [TJ The syst em TV Col was the fir s t syst em to show 
this signal, and iBonnet-Bidaud et al.l (|1985D suggested 
that the periods were consistent with what would be ex- 
pected for a disk that was tilted out of the orbital plane 
and freely precessing with a period of ^4 d. iBarrett et alJ 
(1988) expanded on this and suggested what is now the 
accepted model for the origin of negative superhumps: 
the transit of the accretion stream impact point across 
the face of a tilted accretion disk that precesse s in the ret- 
rograde direction (see IW ood ct al. 2000; Woo d fc Burkd 
l2007t IWood et al.ll2009riFoulkes et al.l 120061) . As in the 
stream source for positive superhumps, the modulation 
results because the accretion stream impact point has 
a periodically-varying depth in the potential well of the 
primary star. 

Finding the term "negative period excess" unnecessar- 
ily turgid, in this work we refer to the period deficit e_ 
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Empirically, it is found that for systems showing both 
positive and negative superhump s that e+/e_ ^ 2 
(|Patter soh 1999; Rcttcr et al.'2002). 

We show in Figure [5] a snapshot from a q — 0.40 simula- 
tion that demonstrates the physical origin of negative su- 
perhumps. At orbit 400, the disk particles were tilted 5° 
about the a:-axis and the simulation restarted. The green 
line in the Figure running diagonally though the primary 
indicates the location of the line of nodes; the disk mid- 
plane includes this line, but is below the orbital plane to 
the right of the line, and above the orbital plane to the 
left of the line. The disk particles are again color-coded 
by luminosity, and the brightest particles are shown with 
larger symbols. The ballistic accretion stream can be fol- 
lowed from the LI point to the impact point near the line 
of nodes. The simulation light c urve is derived from th e 
"surface" particles as described in lWood fc Burkd (|2007D . 
The times of maximum of the negative superhump light 
curve occur when accretion stream impact point is deep- 
est in the potential of the primary and on the side of 
the disk facing the observer. A second observer viewing 
the disk from the opposite side would still see negative 
superhumps, but antiphased to those of the first. 

Having introduced a viable model for positive super- 
humps and their evolution, let us now compare the model 
to the Kepler V344 Lyr photometry. 

3. Kepler PHOTOMETRIC OBSERVATIONS 

The primary science mission of the NASA Discovery 
mission Kepler is to discover and characterize terrestrial 
planets in the h abitable zone of Sun- l ike stars us i ng th e 
transit method (jBorucki et al.l 120101: iHaas et al.l 120101 ). 
The spacecraft is in an Earth-trailing orbit, allowing it 
to view its roughly 150,000 target stars continuously for 
the 3.5-yr mission lifetime. The photometer has no shut- 
ter and stares continuously at the target field. Each in- 
tegration lasts 6.54 s. Due to memory and bandwidth 
constraints, only data from the pre-selected target aper- 
tures are kept. Kepler can observe up to 170,000 tar- 
gets using the long-cadence (LC) mode, summing 270 
integrations over 29.4 min, and up to 512 targets in the 
short-cadence (SC) mode, summing 9 integrations for an 
effective exposure time of 58.8 s. 

There are gaps in the Kepler data streams resulting 
from, for example, monthly data downloads using the 
high-gain antenna and quarterly 90° spacecraft rolls, 
as well as unplanned safe-mode and loss of fine point 
events. For further details of the spacecraft commission- 
ing, target tables, data collection and pro c essing, and 
perfor manc e metrics, see IHaas et all (|2010D . IKoch et all 
(l2010l) . and lCaldwell et al.l (I2010D 

Kepler data are provided as quarterly FITS files by the 
Science Operations Center after being proc essed through 
the standard data reduction pipeline ([Jenkins et al.l 
[20T0I) . The raw data are first corrected for bias, smear 
induced by the shutterless readout, and sky background. 
Time series are extracted using simple aperture photom- 
etry (SAP) using an optimal aperture for each star, and 
these "SAP light curves" are what we use in this study. 
The dates and times for the beginning and end of Q2, 



Q3 and Q4 are listed in Tabled] 

The full SAP light curve for Kepler quarters Q2, Q3, 
and Q4 is shown in flux units in Figure |6l In Figure 2 of 
Paper II we show the full SAP light curve in Kp magni- 
tude units. As noted in Paper II and evident in FigurelHl 
the superoutbursts begin as normal DN outbursts. 

The Q2 data begin at BJD 2455002.5098. For sim- 
plicity we will below refer to events as occurring on, for 
example, day 70, which should be interpreted to mean 
BJD 2455070 - that is we take BJD 2455000 to be our 
fiducial time reference. 

In this paper, we focus on the superhump and orbital 
signals present in the data. The outburst behavior of 
these data in the context of constraining the thermal- 
viscous limit cycle is published separately (Paper II). 

To remove the large-amplitude outburst behavior from 
the raw light curve - i.e., to high-pass filter the data - 
we subtracted a boxcar-smoothed copy of the light curve 
from the SAP light curve. The window width was taken 
to be the superhump cycle length (2.2 hr or 135 points). 
To minimize the effects of data gaps, we split the data 
into a separate file anytime we had a data gap of more 
than 1 cycle. This resulted in 10 data chunks. Once 
the data residual light curve was calculated, we again 
recombined the data into a single file. The results for Q2, 
Q3, and Q4 are shown in Figures[71[Sl andO respectively. 

We also calculated the fractional amplitude light curve 
by dividing the raw light curve by the smoothed light 
curve, and subtracting 1.0. However, as expected, the 
amplitudes of the photometric signals in the residual 
light curve are more nearly constant than those in the 
fractional amplitude light curve. This is because the su- 
perhump signals - both positive and negative - have am- 
plitudes determined by physical processes within the disk 
that are not strong functions of the overall disk luminos- 
ity. 

4. THE FOURIER TRANSFORM 

In Figure [10] we show the discrete Fourier transform 
amplitude spectra for the current data set. We took the 
transforms over 2000 frequency points spanning to 70 
cycles per day. Each transform is of a 5 day window 
of the data, and the window was moved roughly 1/2 day 
between subsequent transforms. The color scale indicates 
the logarithm of the residual count light curve amplitude 
in units of counts per cadence. In Figure [TT] we show a 
magnified view including only frequencies 9.5 to 12.5 c/d 
to better bring out the 3 fundamental frequencies in the 
system. 

Figures [TU] and [TT] are rich with information. The pos- 
itive superhumps (P+ = 2.20 hr) dominate the power for 
days --58-80 and ~162-190. In Figure [TT] we see the that 
time evolution of the fundamental oscillation frequency 
is remarkably similar in both superoutbursts. The dy- 
namics behind this are discussed below in §5.2 where the 
0-C diagrams are presented. 

Once the majority of the mass that will accrete during 
the event has done so, the disk transitions back to the low 
state. This occurs roughly 15 d after superhump onset 
for V344 Lyr. During this transition, the disk source of 
the superhump modulation fades with the disk itself, and 
the stream source of the superhump modulation begins 
to dominate. A careful inspection of Figure [10] shows 
that at this time of transition between disk and stream 
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Fig. 6. — The Kepler Q2, Q3, and Q4 light curve of V344 Lyr in flux units. The outbursts are labeled 1 to 19. 



TABLE 1 

Journal of Observations 



Quarter Start'' End'" 

MJD UT MJD UT 



Q2 55002.008 2009 Jun 20 00:11 55090.975 2009 Sep 17 11;26 
Q3 55092.712 2009 Sep 18 17:05 55182.007 2009 Dec 17 00:09 
Q4 55184.868 2009 Dec 19 20:49 55274.714 2010 Mar 19 17:07 

The start MJD and UT dates are the mid-point of the first cadence of the 
SC time series for each quarter. 

^ The end MJD and UT dates are the mid-point of the last cadence of the SC 
time series for each quarter. 



superhumps, there is comparable power in the second 
harmonic (first overtone) as found in the fundamentaL 
The behavior of the hght curve and Fourier transform are 
more clearly displayed in Figure [T^ which shows 2 days 
of the light curve during the transition period, and the 
associated Fourier transforms. In both cases, the "knee" 
in the superoutburst light curve (see Figure |6] occurs just 
past the midpoint of the data sets. Although the second 
harmonic is strong in transition phase, the pulse shape 
of the disk superhump signal is sharply peaked so the 
fundamental remains prominent in the Fourier transform 
(see Figure [T2|) . 

As can clearly be seen in Figure [11] the orbital period 
of 2.10 hr (11.4 c/d) only becomes readily apparent in 
the Q4 data, starting at about day 200, and it dominates 
the Q4 Fourier transforms. Once identified in Q4, the or- 
bital frequency appears to show some power in the week 
before the first superoutburst in Q2, and between days 
'^130 and the second superoutburst in Q3. Note, how- 
ever, that the amplitude of the orbital signal is roughly 
1 order of magnitude smaller than the amplitude of the 
negative superhump signal, and as much as 2 orders of 
magnitude smaller than the amplitude of the positive su- 
perhump signal. In these data, the orbital signal is found 
only when the positive or negative superhump signals are 
weak or absent. We discuss the physical reason for this 
below. 

Finally, we note that we searched the Fourier transform 
of our Kepler short-cadence (SC) data out to the Nyquist 



frequency of 8.496 niHz for any significant high frequency 
power which might for example indicate accretion onto a 
spinning magnetic primary star (i.e., intermediate polar 
or DQ Her behavior). We found no reliable detection of 
higher frequencies in the data, beyond the well-known 
spurious frequencies present in Kepler time series data 
at multiples of the LC frequency (n x 0.566427 niHz — 
48.9393 c d-i iGilliland et aL|[2010 f). For a fuU list of pos- 
sible spurious frequencies in the SC data, see the Kepler 
Data Characteristics Handbook. 



4.1. The Orbital Period 

The orbital period is the most fundamental clock in 
a binary system In the original Q2 data presented by 
iStill et all (|2010D . the only frequencies that were clearly 
present in the data were the 2.20-hr (10.9 c/d) super- 
hump period and the period observed at 2.06-hr (11.7 
c/d). In Paper I we identified this latter signal as the 
orbital period but discussed the possibility that it is a 
negative superhump period. The Q3 data revealed a 
marginal detection of a period of 2.11 hr (11.4 c/d), and 
this period is found to dominate the Q4 data (see Fig- 
ure [T5)) . The average pulse shape for this signal averaged 
over days 200-275 is shown in Figure [T31 We can now 
safely identify this 2.11 hr (11.4 c/d) signal as the sys- 
tem orbital period, which then indicates that the 2.06 hr 
(11.7 c/d) signal is a negative superhump. 

The orbital period was determined using the method 
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Time (BJD - 2455000) 

Fig. 7. — The residuals of the Kepler Q2 V344 Lyr hght curve after subtraction of a smoothed copy of the SAP light curve. Each panel 
shows approximately 10 days of data (~120 cycles). Note that the vertical scaling of the 2 panels showing the high-amplitude positive 
superhumps {P+ = 2.20 hr, = 10.9 c d~^) differ from the rest (days ~52-72), and that there is a small amount of overlap between 
panels. The negative superhump signal dominates in days ~2-40 (P— = 2.06 hr period, = 11.7 c d~^). The positive superhumps 
dominate the power for days ~58 through the end of the quarter. 



of non-linear least squares fitting a function of the form 
= Asin[27r(i-To)/P]. (4) 
The results of the fit are 

P = 0.087904 ±3 X fO"^ d, 
= 2.109696 ±7 X 10"^ hr, 
To = BJD 2455200.2080 ±0.0006, 
A = 7.8±0.1 e" s-\ 

Note that the amplitude is only roughly 25 mmag - an 
order of magnitude or more smaller than the peak am- 
plitudes of the positive and negative superhumps in the 
system. 



That an orbital signal exists indicates that the system 
is not face-on. The source of the orbital signal of a non- 
superhumping CV can be either the variable flux along 
the line of site from a bright spot that is periodically 
shadowed as it sweeps around the back rim of the disk, or 
the so-called reflection effect as the face of the secondary 
star that is illuminated by the UV r adiation of the disk 
rotates in to and out of view (e.g.. IWarne"^ll995aD . In 
Figure [Til we find that the orbital signal is never ob- 
served when the positive superhumps are present, but 
this is not a strong constraint as the positive superhump 
amplitude swamps that of the orbital signal. 

More revealing is the interplay between the orbital sig- 
nal, the negative superhump signal, and the DN out- 
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Fig. 8. — The residuals of the Kepler Q3 V344 Lyr light curve. The positive superhumps dominate the power for days ~162 through 
the end of the quarter. The negative superhumps dominate the signal during days ~100-160. The vertical scaling of the bottom 2 panels 
differs from the rest. The features that appear during days ~105, 138, 150, and 161 are artifacts of the reduction process. 



bursts. In Q2 and Q3, the orbital signal appears only 
when the negative super hump signal is weak or absent. 
This is consistent with the idea that the addition of mate- 
rial from the accretion stream should bring the disk back 
to the orb ital plane roughly on the mass-replac ement 
time scale ()Wood fc Burke|[2007l: IWood et al.|[20Q9l ). The 
strong negative superhump signal early in Q2 indicates 
a tilt of '^5°, sufhcient for the accretion stream to avoid 
interaction with the disk rim for all phases except those 
in which the disk rim is along the line of nodes. As the 
disk tilt declines, however, an increasing fraction of the 
stream material will impact the disk rim and not the in- 
ner disk - in other words, the orbital signal will grow 
at the expense of the negative superhump signal. This 
appears to be consistent with the data in hand and if 



so would suggest that the orbital signal results from the 
bright spot in V344 Lyr, but the result is only speculative 
at present. 

In Figure [TSl we show the 0-C phase diagram for Porb- 
We fit 20 cycles for each point in the Figure, and moved 
the window 10 cycles between fits. The small apparent 
wanderings in phase result from interference from the 
other periods present, and also appear to correlate with 
the outbursts. We show the 2D DFT for days 200 to 275 
in Figure [TBI Here we used a window width of 2 days, and 
shifted the window by 1 /8th of a day between transforms. 
We show amplitude per cadence. The orbital signal ap- 
pears to be increasing in amplitude slightly during Q4, 
perhaps as a result of the buildup of mass in the outer 
disk after several DN outbursts. The large amplitudes 
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Time (BJD - 2455000) 

Fig. 9. — The residuals of the Kepler Q4 V344 Lyr hght curve. The orbital signal (Porb = 2.11 hr, = 11.4 c d~^) dominates following 
the decay of the postive superhump signal at day ~190. The features that appear during days ~195, 212, 228, 246, and 256 are artifacts 
of the reduction process. 



found for the orbital signal in Figure [15] during outbursts 
17 and 19 (starting days ^246. 5 and 266, respectively) 
are spurious, resulting from the higher-frequency signals 
found on the decline from maximum in each case. As 
discussed below, outbursts 17 and 19 both show evidence 
for triggering a negative superhump signal, and the light 
curve for outburst 19 yields a complex Fourier transform 
that shows power at the orbital frequency, the negative 
superhump frequency, and at 12.3 c/d (1.95 hr). 

4.2. Observed Positive Superhumps 

The light curve for V344 Lyr is rich in detail, and in 
particular provides the best data yet for exploring the 
time evolution of positive superhumps. As discussed 
above, the superhumps are first driven to resonance dur- 
ing the DN outburst that precedes the superoutburst as 



the heating wave transitions the outer disk to the high- 
viscosity state allowing the resonance to be driven to am- 
plitudes that can modulate the system luminosity. Close 
inspection of the positive superhumps in Figures [7] and [8] 
shows that in both cases the amplitude of the superhump 
is initially quite small, but grows to saturation [A ~ 0.25 
mag) in roughly 16 cycles. There is a signal evident pre- 
ceding the second superoutburst (days ~156.5 to 161) 
- this is a blend of the orbital signal and a very weak 
negative superhump signal. 

The mean superhump period obtained by averaging 
the results from non-linear least squares fits to the 
disk superhump signal during the two superoutburst 
growth through plateau phases is P+ — 0.091769(3) d = 
2.20245(8)hr. The errors quoted for the last signifi- 
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Time fBJD - 2455000) 

Fig. 10. — The 2D DFT of the Kepler Data from Q2, Q3, and Q4. Shown here is the logarithm of the amplitude versus time and 
frequency. There are three fundamental periods visible, plus higher harmonics of two of these. The positive superhumps (P-f = 2.20 lir, 
i/+ = 10.9 c d~^) dominate the power for days ~58-80 and ^162— 190. The negative superhumps dominate the signal during days ~100-160 
(P_ = 2.06 hr, V- = 11.7 c d~^) dominate early in Q2, and again in Q3 between the superoutbursts. The orbital signal (Porb = 2.11 hr, 
^ot\> = 11-4 c d~^) is apparent during Q4, and once identified, during the week before the second superoutburst. 




Time (BJD - 2455000) 

Fig. 11. — Magnified view of the 2D DFT for quarters Q2, Q3, and Q4. Shown here is the logarithm of the amplitude versus time 
and frequency. The three distinct frequencies are clearly evident. In days ~2-40 and again in days ~100— 160, the P_ = 2.06 hr period 
{v- = 11.7 c d~^) indicative of the negative superhump dominates. The positive superhumps (P+ = 2.20 hr, I'.f = 10.9 c d~^) dominate 
the power for days ~58-80 and ~162-190. The orbital period (Porb = 2.11 hr, u^ih = 11-4 c d~^) is most clearly apparent in the Q4 data, 
starting about day 200. 
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Fig. 12. — Light curves and associated Fourier transforms dur- 
ing the transition from disk to stream superhumps. As the system 
transitions from disk to stream superhumps in the drop to quies- 
cence, the hght curve is strongly double peaked but non-sinusoidal, 
resulting in power at both and 2u+ (first and second harmon- 
ics). 
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Fig. 13. — The Fourier transform obtained using the Q4 residual 
data. The single dominant peak is at 11.38 c/d (2.11 hr). We 
identify this as the orbital period. 
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Fig. 14. — The average pulse shape for the identified orbital pe- 
riod, averaged over days 200-275 in Q4. Also shown for comparison 
is a sine curve. 
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Time (BJED - 55000) 

Fig. 15. — The O-C phase diagram for the identified orbital pe- 
riod Porb = 2.10 hr. For each point 20 cycles were fit, and each 
window is offset from the previous by 10 cycles. The bottom panel 
shows the amplitude of the fits to the residual light curves. The 
phase is steady, showing only small variations resulting from inter- 
ference with other periods that are present, and that correlate with 
times of outburst. The large amplitudes found during outbursts 17 
and 19 result from the higher-frequency sign al initiated on decline 
from maximum in each case (see Figure [T51l . 

cant digit are the formal errors from the fits summed 
in quadrature. The periods drift significantly during a 
superoutburst, however, indicating these formal error es- 
timates should not be taken seriously. Using the periods 
found for the superhumps and orbit, we find a period ex- 
cess of e+ = 4.4%. We plot the result for V344 Lyr with 
the results from the well-dete rmined systems below the 
period gap listed in Table 9 of [Patterson et aP (2005) in 
Figure [T71 The period excess for V344 Lyr is consistent 
with the existing data. 

In Figures [TH] and [TO] we show the time evolution of 
the mean pulse shape for the first and second superout- 
bursts. To create these Figures, we split the data into 
5-day subsets (^^50 cycles), with an overlap of roughly 
2.5 days from one subset to the next. For each subset we 
computed a discreet Fourier transform and then folded 
the data on the period with the most power. 

The evolution of the mean pulse sha pe is similar to 
results published prev iously (e.g. Patte rson et al.l [20031 : 
IKato et al.ll2009ll2010[ ). however the quality of the Kepler 
data is such that we can test vigorously the model that 
has been slowly emerging in the past few years for the 
origin of the superhump light source, the evolution of the 
pulse shape and the physical origin of late superhumps. 

A comparison of the simulation light curve from Fig- 
ure [T] with the early mean pulse shapes shown in Fig- 
ures [18] and [19] reveals a remarkable similarity, all the 
more remarkable given the very approximate nature of 
the artificial viscosity prescription used in the SPH cal- 
culations and the crude way in which the simulation light 
curves are calculated. 

If the comparison between data and model is correct, 
the SPH simulations illuminate the evolution of the pos- 
itive superhumps from the early disk-dominated source 
to the late stream-dominated source. The signal ob- 
served early in the superoutburst is dominated by disk 
superhumps, where the disk at resonance is driven into 
a large-amplitude oscillation, and viscous dissipation in 
the strongly convergent flows that occur once per su- 
perhump cycle yield the characteristic large- amplitude 
superhumps seen in the top panels of Figures [TH] and 
[T9l After '^100 cycles ('^10 d), a significant amount of 
mass has drained from the disk, and in particular from 
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Fig. 16. — The 2D DFT of the post-superoutburst Q4 data. Shown are amplitude per cadence versus time and frequency. The window 
width was 2 days, and the window was shifted 1/8 day between transforms. The orbital signal is present but weak for most of the quarter. 
Power at a higher frequency is found on the decline from maximum in outbursts 17 and 19. 
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Fig. 17. — Measur ed period excess versus Porb data from 
dPatterson et al.ll2005l) (circles) and for V344 Lyr (triangle). The 
system V344 Lyr is consistent with the trend. 

the driving region. The disk continues to oscihate in re- 
sponse to the driving even after it has transitioned back 
to the quiescent state, but the driving is ofF-resonance 
and the periodic viscous dissipation described above is 
much reduced. Thus, we agree with previous authors 
that the late/quiescent super humps that have been ob- 
served resuh from the dissipation in the bright spot as it 
sweeps around the rim of the non-axisymmetric disk. 

To compute 0-C phase diagrams for each superout- 
burst, we fit a 3-cycle sine curve with the mean period 
of 2.196 hr which yields a relatively constant 0-C during 
the plateau phase. The results are shown in Figures [20] 
and HU The top panel shows the residual light curve 
as well as the SAP light curve smoothed with a window 
width of P+ (135 points). The second panel shows the 
0-C phase diagram, and the third panel the amplitude 
of the fit. Also included in this Figure in the fourth 
panel are the periods of the positive superhumps during 
2-day subsets of the residual light curve obtained with 
Fourier transforms. The horizontal bars show the extent 
of each data window. By differencing adjacent periods, 
we calculate the localized rate of period change of the 
superhumps P^. These results are shown in the bottom 
panel. As perhaps might be expected from the similar- 
ity in the evolution of the mean pulse profile during the 
two superoutbursts, the 0-C phase diagrams as well as 
the evolution of the periods and localized rates of period 



I 

(U 

03 

o 

X 
O 



10 ^ 

Ir 
-10 

10 

i 

-10 I 

5 ^ 

^ 

-5 I 

5 i 

^ 

-5 i 



58.0-63.0 — 
2.21 E 



2 i 
i 



^ 

















-++++- 






68.0- 


-73,0 







-2 ^ 



-2 ^ 

2 














-++++- 






60.0- 


-65,0 



















62.0- 


-67,0 











-++++- 






64.0- 













-++++- 






66.0- 


-71,0 



2.20 

















-++++- 






70.0- 


-75,0 

















-++++- 






72.0- 


-77,0 

















-++++- 






74.0- 


-79.0 

















-++++- 






76.0- 


-81,0 

















-++++- 






78.0- 


-83,0 



2.20 



0.5 



1 



1.5 



Phase 



Fig. 18. — The average positive superhump pulse shapes obtained 
during the first superoutburst. The time spans of each 5-day sub- 
set of data and the folding periods (in hr) are indicated in each 
subpanel. Zero phase is set at the primary minimum. The mean 
pulse shape evolves considerably over the superoutburst. Note the 
changing vertical scale in the top-few panels. 
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Fig. 19. — Similar to Figure lTsl but for the second superoutburst. 
What is particularly interesting is that the evolution of the pulse 
shape is remarkably similar to that seen in the first superoutburst, 
suggesting that the dynamical evolution is also similar from one 
outburst to the next. 

change are also similar. Such diagram s can be illuminat- 
ing in the study of sii perhumps, and iKato et aD ()2009l ) 
and IKato et aLl (j201ClD present a comprehensive popula- 
tion analysis of superhumps using this method. 

When the disk is first driven to oscillation in the growth 
and saturation phase, there is maximum mass at large 
radius, and the corresponding superhump period ('^2.25 
hr) is significantly longer than the mean, yielding a posi- 
tive slope in the 0-C diagram. The rate of period change 
estimated from the first 4 days of data for both super- 



outbursts is P- 



+ 



-8 X 10-4 s s- 



Roughly 10 cycles 



(~ 1 d for V344 Lyr) after the mode saturates with max- 
imum amplitude, sufficient mass has drained from the 
outer disk that the superhump period has decreased to 
the mean, and the superhump period continues to de- 
crease out to i? ~ 100 as the precession rate slows as a 
result of the decreasing mean radius of the flexing, non- 




(BJD - 55000) 



Fig. 20. — The residual light curve (black) in units of 10^ e s ^ 
and smoothed SAP light curve (red) in units of 10* e~ 
{top panel), O-C phase diagram {second panel), amplitudes {third 
panel), period P+ {fourth panel) and localized rates of period 
change P+ {bottom panel) for the common superhumps present 
during the first superoutburst. The observed phase variations re- 
flect real variations in the oscillation period as well as the transition 
in the dominant modulation source for the light curve. The Epoch 
was arbitrarily set to zero for the first datum in this subset of the 
time series. 
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Fig. 21. — Similar to Figure [20l but for the second superoutburst. 
This Figure clearly shows that the phase variations have the same 
character as before, but with slightly larger excursions from the 
mean. 

axisymmetric disk. The period at this time is roughly 
2.19 hr for both superoutbursts, and the rate of period 
change between cycles 30 and 70 which includes the early 
plateau phase before the stream signal becomes impor- 
tant is P+ -1.8 X 10"'' s s-^. 

Between cycles '--^110 and 150, the O-C phase diagrams 
in Figures [20] and [21] show phase shifts of ~0.5 cycles. 
This is the result of the continued fading of the disk 
superhump, and the transition to the stream/late su- 
perhump signal. Careful inspection of the top panels 
of Figures [201 and [m near days 68 and 174 in fact shows 
the decreasing amplitude of the disk superhump, and the 
relatively constant amplitude of the stream superhump. 
By cycle ~ 150 (days ^72 and 176), the disk superhump 
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amplitude is negligible, and all that remains is the signal 
from the stream super hump. The smoothed SAP light 
curve shown in the top panel shows that these times cor- 
respond to the return to the quiescent state during which 
the global viscosity is again low. It is also interesting that 
P+ itself appears to be increasing relatively linearly dur- 
ing much of the plateau phase with an average rate of 
P ~10~^ s~^. At present this is not explained by the 
numerical simulations. It may simply be that this result 
reflects the growing relative importance of the stream su- 
per hump signal on the phase of the 3-cycle sine fit. This 
is almost certainly the case during the period peaks found 
at days ~71 and 175, where we find that the sine fits 
are pulled to longer period by the complex and rapidly 
changing waveform (e.g.. Figure [T^ . 

In the quiescent interval before the first subsequent 
outburst the 0-C diagram shows a concave-downward 
shape indicating a negative P+ ~ —2 x 10~^ s s~^. We 
speculate that the behavior of the 0-C curve in response 
to the outburst following the first superoutburst may in- 
dicate that the outburst may effectively expand the ra- 
dius of the disk causing a faster apsidal precession. Un- 
fortunately, there is a gap in the Kepler data that starts 
just after the initial rise of the outburst following the 
second superoutburst. The value of P+ averaged over 
the the last 2 measured bins for both superoutbursts is 
P+ 3 X 10--* s s-i. 

The measured values of P+ for V344 Lyr are consis- 
tent with t h ose re ported in the extensive compilation of 
IKato et all (|2009( ). To make a direct comparison with 
Kato et al., who calculate P+ over the first 200 cycles 
(i.e., plateau phase), we average all the P+ measure- 
ments out to the drop to quiescence, and find an average 
value of —6 x 10~^ s s~^ for the first superoutburst and 
—9 X 10^^ s s^^ for the second. These values for V344 
Lyr are entirely consistent with the Kato et al. results as 
shown in their Figure 8. 

In Still et' aLl (20iC)f ) we noted that V344 Lyr was un- 
usual (but not unique) in that superhumps persist into 
quiescence and through the following outburst in Q2. 
Other systems that have been observed to show (late) 
superhumps into quiescence more typic ally have short 
orbital periods, including VI 159 Ori ( Patterso n et al.l 
[r995l. ER, UMa (IGao et a l. |[l99 9: Zhao et al. 200^), WZ 
Sge ilPatterson et al.l l2002b|)7an d the WZ Sge-hke star 
V466 And (jChochol et allboiOO . among others. The 
identification of late s uperhumps is a matter of con- 
tention in some cases ()Kato et al.ll2009() . and the post- 
superoutburst coverage of targets is more sparse than 
the coverage during superoutbursts. Thus it is difficult 
to know if post-superoutburst superhumps are common 
or rare at this time. 

4.3. Observed Negative Superhumps 

As noted above in §2.2, the 2.06-hr (11.4 c/d) signal 
that dominates the light curve for the first ^35 days of 
Q2 is now understood to be the result of a negative super- 
hump. This yields a value for the period deficit (Equa- 
tion of e_ = 2.5%. The maximum amplitude at qui- 
escence is A 0.8 mag. Figure shows 10 cycles of the 
negative superhump signal during this time. The inset 
shows the mean pulse shape averaged over days 5 to 25 
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Fig. 22. — The Figure shows 10 cycles of the negative superhump 
signal. The inset shows the average pulse shape over days 5 to 25. 
The signal is sawtoothed with a rise time roughly twice the fall 
time, and appears consistent with the ray-traced SPH simulation 
pulse shapes published by Wood & Burke (2007). 

(roughly 230 cycles). The signal is approximately saw- 
toothed with a rise time roughly twice the fall time. It 
appea rs consistent with the pulse shapes iWood fc Burkd 
()2007[) obtained using ray-trace techniques on 3D simu- 
lations of tilted disks (their Figure 3). Negative super- 
humps dominate the power in days ~2-35 and again in 
days -100-160. 

The signal observed near the beginning of Q2 reveals 
a remarkably large rate of period change - large enough 
that it can be seen in the harmonics of the Fourier trans- 
form shown in Figure [10] as a negative slope towards 
lower frequency with time. A nonlinear least squares 
fit to the fundamental period measured during days 2.5- 
7.5 yields P_ = 2.05006 ± 0.00005 hr. A fit to the data 
from days 22-26, however, yields P_ = 2.06273±0.00005 
hr. The formal errors from non-linear least squares fits 
underestim ate the true errors by as much a n order of 
magnitude (jMontgomerv &: Odonoghiii 119990 . but even 
if this is the case, these two results differ by ^25a. 
Taken at face value, they yield a rate of period change 
ofP_ — 3 X 10^^ s s^^. Similarly, we fit the negative 
superhump periods in two 4-day windows centered on 
days 112.0 and 121.0. The periods obtained from non- 
linear least squares are P_ — 2.0530 ± 0.0002 hr and 
P_ 2.066038 ± 0.00008 hr, respectively, which yields 
P_ — 6 X 10^^ s s~^ over this time span. In their re- 
cent comprehensive analysis of the evolutio n of C Vs as 
revealed by their donor stars, Knig ge^t al.l (|2011f) esti- 
mate that for systems with Poib ~ 2 hr the rate of orbital 
period change should be Porb — — 7 x lO^'' s s~^ (see their 
Figure 11). Clearly the '^2. 06-hr signal cannot be orbital 
in origin. In some negatively superhumping systems with 
high inclinations, the precessing tilted disk can m odu- 
late the mean brightness (e.g. lStanishev et al.ll2002} ). We 
found no significant signal in the Fourier transform at the 
precession period of '^3.6 d. 

In Figure US] we show the results of the 0-C analysis 
for the Q2 data. To create the Figure, we fit 5-cycle sine 
curves of period 2.05 hr to the residual light curve, shift- 
ing the data by one cycle between fits. The shape of the 
0-C diagram is concave up until the peak of the first out- 
burst at day ~28 indicating that the period of the signal 
is lengthening during this time span. The magnitude of 
the negative superhump period deficit is inversely related 
to the retrograde precession period of the tilted disk - a 
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Fig. 23. — The SAP light curve (top panel) O-C phase diagram 
(middle panel) and amplitude (bottom panel) of the negative su- 
perhump signal during days 2.5 to 45.0. A 5-cycle sinusoid was 
fit to the data, moving the v^indow 1 cycle between fits. The am- 
plitude of the signal is large and relatively constant initially, and 
the phase variations are concave upward, indicating a lengthening 
period. During outburst 1, there is a cusp in the O-C diagram, 
indicating a jump to shorter period, and after outburst 1, the am- 
plitude begins to decline. 

shorter precession period yields a larger period deficit. A 
disk that was not precessing at all would show a nega- 
tive superhump period equal to the orbital period. The 
observation that the negative superhump period in V344 
Lyr is lengthening during days ~2 to 27 indicates that 
the precession period of the tilted disk is increasing (i.e., 
the rate of precession is decreasing) . Coincident with the 
first DN outburst (outburst 1) in Q2, there is a cusp in 
the O-C diagram, indicating a jump to shorter period 
(faster retrograde precession rate). The amplitude of the 
signal begins to decline significantly following outburst 
1, and the signal is effectively quenched by outburst 2. 
Note that between days '^28 and 35 the O-C diagram 
is again concave up, although with less curvature than 
before outburst 1. 

We show the 2D DFT of the pre-superoutburst Q2 data 
in Figure [24l Here we used a window width of 2 days that 
was shifted 1/8 day between transforms. We plot the am- 
plitude in counts per cadence. It is evident that outburst 
1 shifts the oscillation frequency, as well quenching the 
amplitude of the signal. Outburst 2 triggers a short-lived 
signal with a period of roughly 11.9 c/d (2.02 hr), and 
outburst 3 appears to generate signals near the frequen- 
cies of the negative and positive superhumps that rapidly 
evolve to higher and lower frequencies, respectively, only 
to fade to to noise background by the end of the outburst. 
Outburst 3 has a somewhat slower rise to maximum than 
most of the outbursts in the time series and is the last 
outburst before the first superout burst, but is otherwise 
unremarkable. This is the only time we see this behavior 
in the 3 quarters of data we present, so it is unclear what 
the underlying physical mechanism is. 

Although much of the Q3 light curve is dominated by 
the negative superhump signal, the amplitude is much 
lower than early in Q2, and in addition there is contam- 
ination from the orbital and positive superhump signals. 
In Figure [25] we show the 2D DFT for the Q3 data be- 
tween days 93 and 162, again showing the amplitude in 
counts per cadence versus time and frequency. We used a 
window width of 2 days that was shifted 1 /8 day between 
transforms. 

In Figure [26l we show the O-C phase diagram obtained 



by fitting a 5-cycle sine curve of period 2.06 hr to data 
spanning days 93.2 to 140.0. The amplitude during this 
time is considerably smaller than was the case for the Q2 
negative superhumps. Before day 106, there appears to 
be contamination from periodicities near the superhump 
frequency of 10.9 c/d that are evident in Figure E5l and 
after day 126 the signal fades dramatically. It was only 
during days 106.5 to 123.2 that the amplitude of the neg- 
ative superhump signal was large enough, stable enough, 
and uncontaminated to yield a clean O-C phase diagram. 
These data lie between outbursts 8 and 9, and comprise 
the longest quiescent stretch in Q3. It can be seen that 
the O-C curve is again concave upward indicating a pos- 
itive rate of period change as calculated above, and the 
bottom panel indicates that the amplitude of the signal 
is increasing during this time span. 

The retrograde precession rate of a tilted ac- 
cretion disk is a direct function of the effec- 
tive (mass weighted) radius of the disk. Several 
groups have studied the precession properties of tilted 
disks ( fPapaloizou & Terquem 1995; Larwood ct al. 199^; 
IJParwood 1998: Tarwood fc Papaloizoui 119971 : lLail il9991. 
Paoaloizou ct aL| (I1997D derived the following expression 
for the induced precession frequency uip of a tilted accre- 
tion disk, 



3 GM2 J dr 
'4 a3 / Y^VLr^dr 



cos (5 



(5) 



where Up is the leading-order term of the induced pre- 
cession frequency for a differentially rotating fluid disk, 
calculated using linear perturbation theory, S(r) is the 
axisymmetric surface density profile and ^{r) the unper- 
turbed Keplerian angular velocity profile, a is the orbital 
separation, M2 is the mass of the secondary, and 5 is the 
tilt of the disk with respect to the orbital plane. The in- 
tegrals are to be taken between the inner and outer radii 
of the disk. 

In a later study of the precession of tilted accretion 
disks. iLarwoodl (|1997l . and see Larwood (1998)) derived 
the expression for the precession frequency of a disk with 
constant surface density S and polytropic equation of 
state with ratio of specific heats equal to 5/3: 



Up 3 f Rq\ 



(6) 



where here VIq is the Keplerian angular velocity of the 
outer disk of radius i?o, and q is the mass ratio. 

The physical interpretation of Equations [5] and [S] is 
that tilted accretion disks weighted to larger radii will 
have higher precession frequencies than those weighted 
to smaller radii. For example, if we have 2 disks with 
the same nominal tilt and total mass, where one has a 
constant surface density and the other with a surface 
density that increases with radius, the second disk will 
have a higher precession rate, and would yield a negative 
superhump frequency higher than the first. A third disk 
with most of its mass concentrated at small radius would 
have a lower precession frequency and yield a negative 
superhump signal nearest the orbital signal. 

In this picture the increasing precession period indi- 
cated by the positive rate of period change for the neg- 
ative superhump signal P_ might at first seem counter- 
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Time (BJD - 2455000) 

Fig. 24. — The 2D DFT of the pre-superoutburst Q2 data. Shown are amphtude per cadence versus time and frequency. The window 
width was 2 days, and the window was shifted 1/8 day between transforms. Most obvious is that the first outburst shifts the power to 
higher frequency, but careful inspection also shows that the frequency evolves towards longer period up to day ~25 (e.g., note the slope of 
the sidelobe at ^ll c/d with respect to the grid line). Outbursts 2 and 3 both appear to generate power on the decline from maximum. 



1.0 

0.5 



V 





12 










>. 




u 


11 


c 


(LI 








w 




\ 


10 


U- 





Too 



no 



Time (BJD - 2455000)' 



150 



IGO 



3200 
2800 
2400 

bnoo 
leno 

1200 

800 

JOO 

to 



Fig. 25. — The 2D DFT of the inter-superoutburst Q3 data. Shown are amplitude per cadence versus time and frequency. The window 
width was 2 days, and the window was shifted 1/8 day between transforms. The negative superhump signal is present until day ~140. 
Higher-frequency power is generated on the decline from maximum of the last 2 outbursts. 

lution of the mass distribution in tilted accretion disks. 

The cause of disk tihs in CVs is still not satisfac- 
torily explained. In the low-mass x-ray binaries it is 
believed that radiation pressure can provide the force 
necessary to tilt the disk out of the orbital plane 
(lPetter son"1977^ 'Iping fc Petter sonI HOOOl: iFoulkes et all 
2006t llvanov & Papaloizou 2008), howeve r this mecha- 



intuitive since the disk is gaining mass at quiescence. 
However, the key fact is that tilted disks accrete most 
of their mass at small radii, since the accretion stream 
impacts the face of th e tilted disk along t he line of nodes 
(|Wood k Burkell2007t IWood et al.l[200g i. The accretion 
stream impacts the rim of the disk only twice per or- 
bit (refer back to Figure [S]). Thus, the effective (mass 
weighted) radius of an accreting tilted disk decreases with 
time, causing a slowing in the retrograde precession rate 
ujp, and an increase in the period of the negative super- 
hump signal P- . 

A detailed analysis of the data, theory, and numerical 
model results should allow us to probe the time evolution 
of the mass distribution in disks undergoing negative su- 
perhumps, and hence the low-state viscosity mechanism. 
The unprecedented quality and quantity of the Kepler 
time series data suggests that V344 Lyr and perhaps 
other Kepler -field CVs that display negative superhumps 
may significantly advance our understanding of the evo- 



nism is not effective in the CV scenario. iBarrett et al.l 
(1_9M) suggested in their work on TV Col that mag- 
netic fields near the LI region might deflect the ac- 
cre tion stream out of the orbital plane, but as noted 
in iWood fc Burkd ()2007f ) the orbit-averaged angular 
momentum vector of a deflected stream would still 
be parallel to the o rbital angular momentum variable. 
iMurrav et al.l ()2002[ ) demonstrated numerically that a 
disk tilt could be generated by instantaneously turning 
on a magnetic field on the secondary star. Although 
their tilt decayed with time (the orbit-averaged angular 
momentum argument again), their results suggest that 
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Fig. 26. — The SAP light curve (top panel) O-C phase diagram 
(middle panel) and amplitude (bottom panel) of the nega tive su- 
perhump signal during days 106.5 to 123.2. As in Figure [23] the 
shape of the O-C curve is concave upward, indicating a positive 
rate of period change. 

changing magnetic field geometries could generate disk 
tilt. A ssuming that the disk v i scosity is con trolled by the 
MRI (Balbus & Hasdey''1998; Balbus 2003), it is plausi- 
ble that differentially-rotating plasmas may also be sub- 
ject to magnetic reconnection events (flares) which are 
asymmetrical with respect to the disk plane, or that dur- 
ing an outburst the intensified disk field may couple to 
the tilted dipole field on the primary star (e.g., Lai 19990 
or the field of the secondary star (Murray et al. 2002). 

With these ideas in mind, the behavior of V344 Lyr 
during outbursts 2, 10, 11, 17, and 19 is tantalizing. 
First, again consider the 2D DFTs from Q2, Q3, and 
Q4 shown in Figures [24 l [25] and [T6 l respectively. In each 
of these cases, there is power generated at a frequency 
consistent with the negative superhump frequency on the 
decline from maximum light. Outbursts 2 and 10 ap- 
pear to excite a frequency of roughly 12 c/d ('^2 hr), 
outburst 17 excites the negative superhump frequency 
for ^3 days, and outbursts 11 and 19 appear to excite 
power at the negative superhump frequency that rapidly 
evolves to shorter frequencies. We show the SAP light 
curves for these outbursts as well as the residual light 
curves in Figure [27l The residual light curves for these 5 
outbursts all appear to show the excitation of a frequency 
near or slightly greater than the negative superhump fre- 
quency that dominates early in Q2. This is about 1/3 
of the normal outbursts in the 3 quarters of Kepler data 
- the other 12 outbursts do not show evidence for hav- 
ing excited new frequencies. Thus, while additional data 
are clearly required and our conclusions are speculative, 
we suggest that these results support a model in which 
the disk tilt is generated by the transitory (impulsive) 
coupling between an intensified disk magnetic field and 
the field of the primary or secondary star. The fact that 
these 5 outburst events yield frequencies near 12 c/d ap- 
pears to support the model that it is the mass in the 
outer disk that is initially tilted out of the plane. 

5. CONCLUSIONS 

We present the results of the analysis of 3 quarters 
of Kepler time series photometric data from the system 
V344 Lyr. Our major findings are: 



1. The orbital, negative superhump, and positive su 
perhump periods are Porb = 2.11 hr, P_ = 2.06 




Time (HJI) - 2455000) 

Fig. 27. — Dwarf nova outbursts 2, 10, 11, 17, and 19, each show- 
ing oscillations on the decline from maximum. The black points 
show the SAP light curves, summed two points for clarity. The red 
points show the residual light curves, also summed two points. 



hr, and P+ = 2.20 hr, giving a positive superhump 
period excess of e+ — AA%, and a negative super- 
hump period deficit of e_ = 2.5%. 

2. The quality of the Kepler data is such that we can 
constrain significantly the models for accretion disk 
dynamics that have been proposed in the past sev- 
eral years. 

3. The evolution of the pulse shapes and phases of the 
positive superhump residual light curve provides 
convincing evidence in support of the two-source 
model for positive superhumps. Early in the super- 
outburst, viscous dissipation in the strongly con- 
vergent flows of the flexing disk provide the modu- 
lation observed at the superhump frequency. Once 
the system has returned to quiescence, the modula- 
tion is caused by the periodically-variable dissipa- 
tion at the bright spot as it sweeps around the rim 
of the still non-axisymmetric, flexing disk. During 
the transition the O-C phase diagram shows a shift 
of ^ 0.5 in phase. 

4. Superoutbursts begin as normal DN outbursts. 
The rise to superoutburst is largely explained by 
the thermal-viscous limit cycle model discussed in 
Paper II. Beyond this luminosity source which does 
a reasonable job of matching the lower envelope of 
the superoutburst light curve, there is additional 
periodic dissipation that generates the superhump 
signals. The sources of the periodic dissipation are 

(i) the strongly convergent flows that are generated 
once per superhump cycle as the disk is compressed 
in the radial direction opposite the secondary, and 

(ii) the variable depth of the bright spot as it 
sweeps around the rim of the non-axisymmetric os- 
cillating disk. 

5. Numerical experiments that individually isolate the 
two proposed physical sources of the positive super- 
hump signal yield results that are broadly consis- 
tent with the signals in the data. 
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6. The positive superhumps show significant changes 
in period that occur in both superoutbursts. The 
average P+ ~ 6 x 10^^ s s~^ for the first super- 
outburst and /+ ~ 9 x 10~^ s s~^ for the second 
are consistent with hterature results. The data re- 
veal that P+ itself appears to be increasing rela- 
tively linearly during much of the plateau phase 
at an average rate for the two superoutbursts of 
P ~10-9 s-i. 

7. The negative superhumps show significant changes 
in period with time, resulting from the changing 
mass distribution (moment of inertia) of the tilted 
disk. As the mass of the inner disk increases be- 
fore outburst 1, the retrograde precession period 
increases, consistent with theoretical predictions. 
These data are rich with unmincd information. 

8. Negative superhumps appear to be excited as a di- 
rect result of some of the dwarf nova outbursts. We 
speculate that the MRI-intensified disk field can 
couple to the field of the primary or secondary star 
and provide an impulse that tilts the disk out of 
the orbital plane. Continued monitoring by Kepler 
promises to shed light on this important unsolved 
problem. 

The system V344 Lyr continues to be monitored at 
short cadence by the Kepler mission. It will undoubt- 
edly become the touchstone system against which obser- 



vations of all other SU UMa CVs will be compared, as 
the quantity and quality of the time series data are un- 
precedented in the history of the study of cataclysmic 
variables. The Kepler data for V344 Lyr promise to re- 
veal details of the micro- and macrophysics of stellar ac- 
cretion disks that would be impossible to obtain from 
ground-based observations. 
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